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Abstract 

We study a model of asymptotically free theories with bound states using the similarity 
renormalization group for hamiltonians. We find that the renormalized effective hamiltonians 
can be approximated in a large range of widths by introducing similarity factors and the 
running coupling constant. This approximation loses accuracy for the small widths on the 
order of the bound state energy and it is improved by using the expansion in powers of 
the running coupling constant. The coupling constant for small widths is order 1. The 
small width effective hamiltonian is projected on a small subset of the effective basis states. 
The resulting small matrix is diagonalized and the exact bound state energy is obtained 
with accuracy of the order of 10% using the first three terms in the expansion. We briefly 
describe options for improving the accuracy. 

PACS Numbers: ll.10.Gh 
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1. INTRODUCTION 



So far, we do not have a precise theoretical description of the bound states of quarks and 
gluons in QCD which could simultaneously explain the parton model and the constituent 
quark model of hadronic structure. In particular, QCD is asymptotically free and the per- 
turbative running coupling grows at small momentum transfers beyond limits. This rise 
invalidates the usual perturbative expansion in the region of scales where the bound states 
are formed. 

Ref. suggested a light-front hamiltonian approach to this problem which is based on 
the calculation of the effective hamiltonians using the similarity renormalization group. || 
H Besides the calculation of the hadronic spectrum, one of our purposes is to obtain the 
quark bound state wave functions that can be used in the parton models of large momentum 
transfer processes. An alternative approach is the lattice gauge theory which is making 
progress but does not easily yield such wave functions. M Other recent research in the 
renormalized light-front hamiltonian approach to QCD is described in Ref. |J. 

Wegner |J proposed a flow equation for hamiltonians in solid state physics which is of 
the same kind as in the similarity renormalization group. Wegner's equation is based on 
an explicit form for the generator of the similarity transformation and corresponds to the 
Gaussian similarity factor with a uniform width. 

This paper describes a numerical study of the key elements of the hamiltonian approach 
in a simple matrix model which is asymptotically free, contains a bound state and can be 
diagonalized exactly. We check the accuracy of different approximations in perturbation 
theory by comparison with the exact solution. The exact solution for the renormalization 
group flow of the effective hamiltonians in the matrix model is obtained numerically using 
Wegner's equation. 

We start from the general assumption that the asymptotically free theories have many 
degrees of freedom which are characterized by different scales of energy as measured by cer- 
tain Hq. Then, we define matrix elements of the initial interaction hamiltonian Hj between 
eigenstates of H . 

The model we study can be alternatively derived by discretization of the 2-dimensional 
Schrodinger equation with a potential of the form a coupling constant times a 5-function. []7| 
The continuous version of the 2-dimensional model has been studied by many authors. The 
discretized version was analyzed in Ref. || using the exact solution to Wegner's equation. 
It was shown there that the Wegner equation has the renormalization group interpretation. 

This paper is organized as follows. In Section 2 we describe the model. The parame- 
ters are chosen in a way which will make it clear that the method of solution we use may 
have a wide range of other applications. Numerical results for the effective hamiltonians 
are presented in Section 3. We discuss the approximation based on the similarity factors 
and the running coupling constant and we describe results obtained in perturbation theory. 
Section 4 concludes the paper with a discussion of some options for improvements in pertur- 
bative calculations of the effective hamiltonians in the range of widths near the bound state 
formation scale. 
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2. MODEL 



The hamiltonian H is assumed to have a finite discrete set of eigenstates, 

H \i > = Ei\i > . (2.1) 
The eigenstates are orthogonal and normalized, 

< i\j > = 5 ij . (2.2) 

The dynamics of states in the space spanned by this set is defined by the interaction hamil- 
tonian, Hj, whose matrix elements are assumed to be 

<i\Hj\j> = -g^E~E~, (2.3) 

where g is a dimensionless coupling constant. The whole hamiltonian is denoted by H, 
H = H + Hj. 

In the current study, we choose the eigenvalues of Hq in the form 

Ei = b 2i . (2.4) 

In the numerical calculations we use b = y/2. The integer power i ranges from M to N. 
The integer M is considered to be large and negative. The lower, infrared bound on the free 
energies is 2~ M . The integer iV is considered to be large and positive. The upper, ultraviolet 
free energy bound is given by 2 N . In the numerical study, we use M = —21 and N = 16. 

For the purpose of analogy to QCD, we adopt the convention that the energy equal 1 
corresponds to 1 GeV, although the units of energy are arbitrary. Thus, the ultraviolet cutoff 
corresponds to 65 TeV and the infrared cutoff corresponds to 0.5 eV. 

With the above choices, the hamiltonian H is a 38 x 38 matrix. For the coupling con- 
stant g > 1/38, the matrix has one negative eigenvalue and 37 positive eigenvalues. The 
coupling constant is adjusted to obtain the negative eigenvalue equal —1.00000000 or, in our 
convention, —1 GeV. Namely, g = 0.06060600631. The many digits are given for readers 
interested in the numbers. 

The eigenstate with the negative eigenvalue corresponds to the s-wave bound state in the 
continuum Schrodinger equation with the 5-potential. The positive eigenvalue eigenstates 
correspond to the s-wave scattering states. We refer the reader to the work of Jackiw for 
details. The only new steps required are the replacement of the continuous energy scale in 
Ref. |7| by the discrete one in Eq. (2.4) and the introduction of the infrared and ultraviolet 
cutoffs. These steps are described in Ref. ||. 

To verify the ultraviolet renormalizability of the model and its infrared convergence, we 
have studied a set of exact results for different numbers M and N. The same qualitative 
results as for M = —21 and N = 16 can be obtained already for — M = N = 8. We 
have checked the renormalizability by varying the ultraviolet limit N from 4 to 16 and we 
verified the infrared convergence by varying the limit M from —4 to —21. We have studied 
in detail the case — M = N = 12. The eigenvalues in the latter case are almost the same as 
the corresponding eigenvalues in the case M = —21 and iV = 16 with the accuracy better 
than 1% for the extreme eigenvalues. The intermediate eigenvalues match with much higher 
accuracy, correspondingly. For example, eigenvalues order 1 have the same 5 significant 
digits in both cases. 
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We calculate effective hamiltonians using the similarity renormalization group equations 
in the differential form. The effective hamiltonians are parametrized by their width in energy. 
The width is denoted by A. 

We use Wegner's flow equation || which provides a very elegant definition of the similarity 
transformation with a uniform band width and an explicit expression for the generating 
matrix. The effective hamiltonian matrices, 7i = TC(s), are parametrized by the parameter 
s. s ranges from to oo. Wegner's original notation for this parameter was i instead of s. 
It will be shown that the hamiltonian width is given by A = 1 / ^/s. 

The effective hamiltonian is divided into two parts, 

U = V + H-V = V + V. (2.5) 
T> is the diagonal part of the effective hamiltonian matrix, 

where V m = H mm . V is the effective interaction. It is equal to the off-diagonal part of the 
matrix 7i. The Wegner flow equation is || 

^ = [[V,H),H], (2.7) 

with the initial condition 

H(0) = H . (2.8) 

Equations (2.7) and (2.8) ensure that Ti,(s) is a unitary transform of the initial Hamiltonian 
H. In terms of the matrix elements, we have 

= -(V m -V n ) 2 V mn + J2( V m + V n -2V e )V me V in , (2.9) 

as g 

Equation (2.9) can be approximately solved for a small coupling constant g by keeping 
only terms order 1 and g. Namely, 



ftmn(s) = E m 5 mn - g\JE m E n exp[-s(£? m - E n ) J . (2.10) 

In this approximation, T> m = (1 — g)E m . It is clear that the parameter s and the width A 
of the similarity renormalization scheme are simply related, s = A -2 . The similarity factor 
is a Gaussian function. 

It is well known that the coupling constant g must depend on the upper energy bound, 
2 N , in order to eliminate the logarithmic divergences in higher orders when iV — > oo. Ref. 
0] demonstrated that the Wegner flow equation has the standard renormalization group 
interpretation. Lowering the width A = s -1 / 2 in the discrete model is similar to lowering 
the upper bound on the energies of the interacting states. Therefore, in the higher order 
analysis, we replace the coupling constant g in Eq. (2.10) by the running coupling g(s) and 
use the expansion in g(s) to remove the logarithmic divergences. We will gradually change 
the notation from the mathematical parameter s to the hamiltonian energy width A in our 
parametrization of the effective hamiltonians, including the running coupling constant. 
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The small coupling and large width approximation for g(s) can be derived in the following 
way. E m and E n are very small numbers for m and n close to M. Therefore, Eq. (2.9) for 
m and n close to M, g <C 1 and s close to 0, reduces to 

f = (2.11) 

Integration of Eq. (2.11) for our choice of the model parameters gives the approximate 
running coupling constant 

UX) = l-0(17.4-logA71og4) • (2 - 12) 

A is taken in units of GeV. The number 0.4 results from Eq. (2.11), because the contributions 
of terms with i close to log A 2 / log 4 are smaller than 1. 

Eq. (2.12) implies that the family of the effective theories exhibits typical asymptotic 
freedom behavior: the coupling gets smaller when the effective cutoff A is large. 

The effective hamiltonians can be written now as 



H mn (A) = E m 5 mn - g a (X) \jE m E n exp \-[E m - E n ] /A J + corrections . (2.13) 

Here, V m {\) = [1 — g a (X)]E m + corrections. 

Eq. (2.13) demonstrates the utility of Wegner's equation. However, the uniform width 
of the effective hamiltonian distinguishes this solution from the widening band structure 
in the similarity scheme from Refs. and [^]. The widening of the band is useful in 
high order perturbation theory. A whole class of generalized Wegner equations for the 
hamiltonian matrix elements can be written which allow the high energy widening of the 
effective hamiltonian width. Namely, 

H = ^ imh n x ) . (2.i4) 

In the original Wegner case, E{7i} = [D,?i\. We have made calculations in the model 
using different formulae for the operation F. For example, we used [F{Ti]] mn = 9[\A mn \ — 
x](\A mn \ - x) k A mn Hmn, where A mn = (V m - V n )/(V m + V n ) and a; is a function of A such 
that 1 > x(A) > x > 0, (see also Ref. [|K]]). Nevertheless, we limit the present paper 
discussion to the application of the original Wegner equation, for simplicity and because no 
clear advantage of the more general equations over the Wegner equation has been visible in 
the numerical studies we performed so far. 

The structure of Eq. (2.10) demonstrates that the approximate solution including the 
running coupling cannot be obtained in the first order perturbation theory. Eq. (2.11) shows 
that the running coupling constant approximation can be obtained in perturbation theory 
if one keeps terms order g 2 . The approximate running coupling ^ a (A) is given in Eq. (2.12). 

The question is how large are the corrections indicated in Eq. (2.13). The next Section 
describes our numerical study of the model. 
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3. NUMERICAL RESULTS 



The goal of our numerical study is to understand the structure and determine the useful- 
ness of the effective hamiltonians for the calculation of the bound state energy. We want to 
find out if perturbatively calculated effective hamiltonians of small width (of order 1 GeV) 
can reproduce the bound state eigenvalue. We are interested in perturbation theory because 
it is the tool we can use in the hamiltonian approach to QCD. We aim to make the width A 
as small as possible because the smaller is the width the smaller basis of effective states is 
required to calculate the bound state eigenvalue. 

The questions how small can be the width of the perturbatively calculated effective 
hamiltonians and how well the hamiltonians can reproduce the bound state properties, are of 
principal interest because the perturbative running coupling constant grows in asymptotically 
free theories when the width gets smaller and at some unknown point perturbation theory 
becomes useless. 

The approximate result for the running coupling constant g a (X) in Eq. (2.13) can be 
compared to the exact value of the running coupling. The exact solution is obtained by 
the numerical integration of the flow equation, i.e. Eq. (2.7). The numerical computations 
were performed using the Runge-Kutta integration algorithm of rank 4. || The results 
were obtained using two independent alghorithms. The hamiltonians were cross-checked 
using various theoretical conditions such as the hermiticity and width-independence of the 
eigenvalues. 

The exact running coupling, g(X), is defined as 

g(X) = - n M ,M+i(X)/y/E M E M+1 . (3.1) 

Fig. 1. shows the running coupling constants g and g a as functions of the width A. It is 
visible that the approximate solution blows up in the flow before the effective hamiltonian 
width is reduced to the scale where the bound state is formed. 

The bound state formation scale is defined using the relevant effective hamiltonian matrix 
element fx which is defined as 

/2(A) = W-i_i(A) - 0.5 GeV . (3.2) 

/2 becomes equal —1.5 GeV when the diagonal matrix element Ti-i -i of the effective hamil- 
tonian becomes equal to the bound state eigenvalue, —1 GeV. The width A where the bound 
state eigenvalue appears on the diagonal is equal about 0.5 GeV. /2 stays constant for smaller 
A. This width scale (order 1 GeV) is called the bound state formation scale. The exact result 
of integrating Eq. (2.7) gives the function /2(A) which is plotted in Fig. 1. 

Fig. 1 also shows that the exact effective coupling constant is close to the one in Eq. 
(2.13). This is a hint for that the corrections indicated in Eq. (2.13) are small for A larger 
than order 16 GeV. The effective hamiltonians for smaller A may still have similar structure 
but the running coupling is not given by Eq. (2.12). 

The most important feature visible in Fig. 1 is that the exact effective coupling constant 
does not grow unlimitedly. This encourages us to use expansion in the running coupling 
constant in calculations of the effective hamiltonians. 

On the other hand, the diagonal matrix elements T> m = (1 — g)E m with small m become 
negative when g grows above 1. |§ Thus, the absolute energy order of states is reversed. 
The diagonal matrix elements for states corresponding to lowest E m become negative but 
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they are small in modulus. At the same time, the diagonal matrix elements for states 
corresponding to larger E m become more negative. Therefore, when g grows above 1, the 
states corresponding to some originally intermediate energy range become much lower on 
the energy scale than all other states. 

This inversion feature reverses the role of the bilinear terms. The coupling stops growing 
and it begins to drop. At the same time the bound state eigenvalue is localized in the diagonal 
matrix element 7i-i-\. This matrix element corresponds to the state which appears to be 
at the bottom of the energy scale. The bound state dynamics is located in matrix elements 
at the center of the effective hamiltonian matrix instead of the lowest indices corner, as one 
might expect if the lowest momentum scales were important in the bound state formation. 
Further renormalization group flow reduces g below 1 and successively establishes small 
positive eigenvalues on the diagonal. In other words, the negative bound state eigenvalue 
is settled before many low positive eigenvalues are. The bound state eigenvalue appears 
almost independently of what happens at the bottom of the positive spectrum. This can be 
verified by changing M. The hamiltonian width A limits couplings of a given state with other 
states to a limited number of states whose diagonal matrix elements are within the range 
order A around the given state diagonal matrix element. For b 2 = 2 only a small number of 
states participate in the interaction of states with the diagonal matrix elements larger than 
A. Therefore, when A drops below the bound state formation scale, the small energy (small 
diagonal matrix element) states are decoupled from the bound state dynamics. 

The rise of the coupling above 1 and the inversion feature suggest that even in the full 
hamiltonian flow the perturbative expansion is useless at the bound state formation scale. 
But the coupling approaches 1 from below rather smoothly. We can ask for how small A the 
effective hamiltonian can be reliably calculated in perturbation theory. 

Scaling symmetry in the model and Fig. 1 suggest a large range of validity of perturbation 
theory down to at least the width order 16 GeV where the formula (2.12) begins to fail. 
Namely, the hamiltonians in Eq. (2.13) possess the discrete scaling symmetry: two effective 
hamiltonians for two values of the width A differing by the factor b~ 2 have matrix elements 
related by the shift % — > % — 1 and j — > j — 1 and replacement of g(X) by g(b~ 2 X). Thus, 
the flow of the effective hamiltonians is reduced to the change of the width and the coupling 
constant for as long as one can neglect the boundary effects due to finite N and M. Since 
the coupling in Eq. (2.12) is obtained keeping terms order g 2 and it runs correctly down to 
about 16 GeV, one can expect that Eq. (2.13) describes the effective hamiltonians correctly 
down to that scale. 

Fig. 2 shows the bound state eigenvalue of approximate effective hamiltonians in ratio 
to the exact value —1 GeV, as a function of the width A. Ratios formed this way will be 
used as measures of the accuracy of effective hamiltonians throught this work. 

Three approximations are shown in Fig. 2. The one denoted by E results from diago- 
nalization of hamiltonians given by Eq. (2.13) with corrections set equal and inserting the 
exact value of g(\) for g a (A). The curve labeled by V is obtained in the same way except 
for one modification that in the exponent E, m and E n are multiplied by the factor 1 — g(\). 
This multiplication replaces E m and E n in the exponent by V m and V n , respectively. The 
intermediate curve marked g is obtained by using T> — (1 — g)E in the exponent, i.e. the 
initial coupling instead of the running one. 

The three curves show strong deviation from the exact eigenvalue at the beginning of the 
flow. This is caused by the boundary effect due to the cutoff N. The curves also deviate 
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from 1 at small width. This deviation shows the limited validity of the approximation. It is 
also clear that the gaussian similarity factor with free energies gives a better approximation 
than the same factor with the diagonal matrix elements. 

Fig. 2 suggests that a few terms of expansion in the running coupling (e.g. only terms 
order 1 and g are used in the case of curve E) may produce effective hamiltonians whose 
bound state eigenvalue accuracy is order 10% and whose width is smaller than 16 GeV. The 
amazing result is that by simply introducing the similarity function and varying the coupling 
one can reduce the width of the hamiltonian by a factor order 4000 and make a small error. 
Other eigenvalues are more accurate with the exception of eigenvalues order or larger than 
A for which the approximation is not expected to work. However, the width is still quite far 
from the bound state formation scale. The key question now is how far down in the width 
one can get using an expansion in the running coupling constant. 

The first comment we wish to make concerning perturbation theory is that the direct 
expansion in powers of the initial coupling g is not useful for calculating effective hamiltonians 
with small width. This is illustrated in Fig. 3 by the plot of the ratio of the effective 
hamiltonian bound state eigenvalue to the exact result for hamiltonians of width A calculated 
using expansion into a series of powers of g up to 1, 2, 3 and 4. Note that the analogy between 
the model and QCD is such that g in the model corresponds to the strength of the two-particle 
interaction which is of the second order in QCD. Therefore, terms order g here correspond 
to terms order gq CD , terms order g 2 correspond to terms order gq CD etc. Fig. 3 clearly 
demonstrates that perturbative expansions in terms of the canonical coupling constant in 
the initial hamiltonian are not suitable for applications in the bound state dynamics. 

The reason we stress this fact here despite that such result might be expected in the 
lagrangian approach, is that we study a hamiltonian approach. The hamiltonian approach is 
different in many respects from the familiar perturbative lagrangian approaches. [ lj In the 
light-front hamiltonian approaches, it is often silently assumed that one can analyse nonper- 
turbative aspects of the light-front QCD hamiltonian dynamics using canonical hamiltonian 
terms and neglecting the running coupling effects. Fig. 3 warns us not to do so. 

In contrast to Fig. 3, Fig. 4 shows results obtained from effective hamiltonians calculated 
in the second order expansion in terms of the effective coupling constant (?(Ao), f° r a se ^ °f 
choices of Ao- The curves are marked by the value of Ao in GeV. For comparison, we include 
the second order result of the expansion in terms of g which is marked by the sign oo and 
equals to the curve marked 2 in Fig. 3. The actual value of g which is used to evaluate the 
hamiltonians 7i(X) is equal to the exact value of g(Xo) in the model. In a theory where an 
exact solution is not known, the exact coupling constant g(Xo) is not known and its value 
must be fitted in TC(X) at a useful value of A. Varying of A should not cause significant 
changes if the approximation and the fit are nearly exact. 

The expansion of the renormalization group flow in powers of an effective coupling g(Xo) 
is done in the following way. One expands the flow in powers of the initial g and computes 
g(X ) as a series in g using Eq. (3.1). The latter series is inverted and g is calculated as 
a series in g(Xo). Then, the whole family of effective hamiltonians parametrized by A is 
calculated in expansion of powers of g(Xo) by substituting the inverted series into the known 
expansion in powers of g. 

Thus, the effective hamiltonians are calculated using expansion in powers of g{Xo). They 
are diagonalized and their bound state eigenvalues are plotted as functions of A, for every 
single value of Aq we choose. It is expected that an expansion in powers of g(Xo) works only 
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for effective hamiltonians with widths A in the vicinity of Ao- This feature is clearly visible 
in Fig. 4 for A = oo, 512, 32, 8, 2 and 1 GeV. The arrows in Fig. 4 show points where 
A = A . 

The next question is what happens in the higher orders of perturbation theory. This 
question can be only partially answered on the basis of the limited numerical studies we 
have performed. Therefore, we limit the discussion here to the analysis of the lowest orders 
in the perturbative expansion for the key matrix element of Eq. (3.2) and the bound state 
eigenvalue. 

Firstly, note that Eq. (2.12) can be viewed as a result of summing a geometric series. It 
is clear that when A approaches 2 GeV, the coefficients in the expansion of g a in powers of 
g grow as powers of N. Therefore, one should expect that when g > l/N the series will not 
converge. The value of g we are using in the model is only a little bit smaller than 1/16, 
as required by the condition for the bound state energy to be —1 GeV. For A ^> 2 GeV the 
coefficients are powers of much smaller numbers than N and one can expect a small number 
of terms to reproduce the running coupling value, but when the width gets small, the number 
of terms in the expansion must become very large to reach the approximate result. When 
one switches to the expansion in the exact running coupling constant, the number of needed 
terms is unknown. We will discuss the lowest 6 terms only. 

Secondly, it follows from Fig. 1 that the exact solution for g(X) is limited and always 
stays below 1.1, in a dramatic distinction from the approximate running coupling result 
which diverges below 4 GeV. This suggests that the perturbative expansion for the full 
hamiltonian may be extendible beyond the barrier around 4 GeV. We show below it does 
happen so for the first three terms in the expansion. 

The difficulty we encounter with terms of higher order than g 2 can be seen still using Eq. 
(2.12). Suppose A = 16 GeV. We see in Fig. 1 that the second order perturbative result for 
g a (16GeV) is quite close to the exact coupling constant ^(16 GeV). So, let us approximate 
the function g(X) by the function g a (\) for A > 16 GeV. Then, using Eq. (2.12) at A = 16 
GeV, one obtains the inverse series 

g = g - (13Ag) 2 + (13Ag) 3 - (13.4 ~gf + ... , (3.3) 

with the alternating sign. The first 6 actual coefficients which replace the successive powers 
of 13.4 in the full perturbative expansion are: 1.0, -13.42, 177.3, -2309, 29752 and -379277. 
For example, the approximation in Eq. (3.3) deviates from the actual coefficients by about 
5% for the coefficient of g 4 and about 20% for g 6 . 

Now, let us calculate the matrix element /2(16GeV) of Eq. (3.2) as a power series in 
g. The actual coefficients are -0.50, -6.71, -88.6, -1153, -14832 and -188759. Inserting 
the expansion from Eq. (3.3) one obtains the coefficients in /t(16GeV) in front of the 
powers of #(16 GeV), from the 1st to 6th, equal -0.5, 0.0008, 0.04, -0.14, 1.2 and -17.3, 
correspondingly. One can easily test the series and see that it provides a good convergent 
approximation to the jx(g) for very small values of g, such as 0.001 or even 0.01. 

The exact value of the running coupling constant, g (16 GeV), is of the order 0.27. For 
this value of ~g one obtains the result /2(16 GeV) = [-0.13500 + 0.000057+0.00079-0.00074 + 
0.00174 — .00671] GeV. This sum does not include terms order 1. The first term comes from 
the order g, second from g 2 etc.. 

It is clear now that for the first three terms order g k with k — 0, 1 and 2 we can expect 
convergence of the effective hamiltonian calculation while the higher orders may destroy it 
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if one keeps only a few terms. 

The need for many orders in the expansion to reach convergence beyond the first three 
terms in the expansion can be suggested already on the basis of Fig. 1 . Namely, we see that 
the perturbative switching of g a (X) from +00 to —00 is smoothed out in the full renormal- 
ization group flow of g(X). However, the turn of the curve g below 1 GeV is sharp and a 
large number of terms is required to reproduce the whole curve. 

It is also visible, that the successive terms may cancel out. For example, the fourth order 
term may remove a large contribution of the third order term. If the numerical cancelation 
appears for many digits, the round off errors require a close examination. 

We have not performed extended studies of the model beyond the lowest order terms. 
Therefore, we limit ourself to the presentation of Figs. 5 and 6 which illustrate what happens 
up to terms order g A . The effective hamiltonian TL{\) is expanded in a series of powers of 
g(X) (i.e. the running coupling at the same value of A). 

Fig. 5 shows the accuracy of the perturbatively calculated effective hamiltonians as 
measured by the matrix element /2(A) from Eq. (3.2) in ratio to the exact value. The 
expansion including terms order 1 and g works with 10% or even better accuracy down to 
the width 1 GeV. It is also visible, that the expansion including terms order g 2 produces 
accuracy on the order of 3% or better down to the width order 2 GeV but the accuracy 
drops down significantly near 1 GeV. The third power term produces a considerable drop in 
the accuracy which is then partially counterbalanced by the inclusion of the fourth power 
term. The oscillation pattern needs further studies to explain. We have not performed such 
studies. 

Fig. 6. shows the accuracy of the bound state eigenvalue obtained from diagonalization 
of the effective hamiltonians 7i(A) expanded in powers of the running coupling g(X) including 
powers up to 1, 2, 3 and 4, respectively. The curve 1 for the first order calculation matches 
curve E in Fig. 2 away from the ultraviolet (left) boundary region. This feature confirms 
the expectation based on Figs. 1 and 2 that a low order perturbative expansion in the run- 
ning coupling constant can accurately reproduce the effective hamiltonians at small width. 
The curve 2 represents the result of the second order expansion and shows a considerable 
improvement. The accuracy is about 10% or better down to the width 1 GeV which is at the 
edge of the bound state formation scale. The second power of g in the model is analogous to 
the fourth power of the coupling constant in QCD. Thus, the analogous calculation in QCD 
requires the fourth order expansion of the effective hamiltonians in powers of the strong 
interaction coupling constant. 

It is visible in Fig. 6 that for the small values of g the perturbative expansion is con- 
vergent. However, as expected, at small widths the powers 3 and 4 of the coupling constant 
appear with large coefficients. The curve including terms up to the third power of g(\) falls 
down quite low when we substitute the exact value of the coupling. 

The fourth order curve labeled 4 is ended at the bottom edge of the figure to avoid 
overlap with curves 1 and 2. In fact, the curve 4 continues down to about 0.53 at A = 16 
GeV, shoots back up to 0.94 at about 12 GeV, deeps down again to about 0.5 at 8 GeV 
and skyrockets right after 8 GeV crossing 1 and reaching about 100 at 1 GeV. This erratic 
behavior is clearly correlated with the pattern visible in Fig. 5. 

Despite the slow convergence problem which requires further study, Figs. 1 to 6 illustrate 
the striking feature of the model that the bound state eigenvalue can be obtained from 
diagonalization of an effective hamiltonian with the width order 1 GeV with 10% accuracy. 
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The remaining question is how small the space of states can be on which one can project 
the narrow effective hamiltonian and reproduce the bound state eigenvalue by diagonalization 
of the projected matrix. The wave function is expected to be reproduced with a similar 
accuracy. 

The answer is provided in Table 1. The table contains the ratio of the bound state 
eigenvalue obtained by the diagonalization of a small window matrix 7Y(A) whose matrix 
elements are the same as the matrix elements of 7i(A) in the small window but they equal 
everywhere outside the window. The indices of the window range from rh to n including 
the limiting values. The bound state eigenvalue of the window is divided by the bound 
state eigenvalue of the whole effective hamiltonian. Table 1 gives the results obtained from 
H(l GeV) which is calculated in second order expansion in the running coupling constant 
g(l GeV). The bound state eigenvalue of the whole effective hamiltonian 7i(l GeV) is equal 
—0.8902 GeV (the exact value is —1 GeV). The entries in Table 1 show that the small window 
matrix easily reproduces this result with a relatively high accuracy. This is a spectacular 
feature of the method and the model. For the coupling constant g(l GeV) is equal 1.05 and 
still the bound state eigenvalue accuracy one obtains from the small window hamiltonians 
calculated in second order perturbation theory is on the order of 10%. 

Table 1 demonstrates that our hamiltonian approach can be used to calculate a 10 x 10 or 
even 5x5 matrix whose lowest eigenvalue reproduces the full theory bound state eigenvalue 
with accuracy order 10 to 20%. This model result can be viewed as encouraging to pursue a 
similar strategy in case of QCD. However, the slow down or absence of convergence beyond 
the second order expansion near the bound state formation scale, require improvements. 



4. CONCLUSION 

We have studied two basic approximations which may be of help in the hamiltonian 
calculations of bound state properties in asymptotically free theories. The first one can be 
briefly described as constructing effective hamiltonians by introducing the similarity fac- 
tors and adjusting coupling constants. We provide the definition of the asymptotically free 
running coupling constant in the hamiltonian approach. The second approximation is an 
expansion into a series of powers of the running coupling constant. In both cases, the bound 
state eigenvalues and eigenstates of the full theory are found by diagonalizing the effective 
hamiltonians. 

Our model study suggests that one may hope to calculate effective hamiltonians down 
to the similarity width which is close to the bound state formation scale. The coupling 
constant growth is limited. Moreover, the effective hamiltonian can be diagonalized in a 
limited subspace of states which dominate the bound state dynamics, instead of using the 
full basis. The small window hamiltonian reproduces the whole effective hamiltonian bound 
state eigenvalue with accuracy order 10% or better. 

The model exhibits an inversion of the energy hierarchy of states when the coupling 
constant becomes slightly larger than 1. On the one hand, this feature is welcomed because 
it decouples the small momentum scales from the bound state dynamics. On the other 
hand, this feature appears at the coupling larger than 1 and it is beyond the reach of a 
simple perturbation theory. One can ask if a similar difficulty must exist in QCD. The 
expected answer is no. For one may hope that gluons effectively obtain masses already 
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at small values of the coupling constant through the gluonic couplings. These masses (or 
potentials) may set the states with gluons appart in energy when the width gets small before 
the coupling becomes too large for being treated in perturbation theory. 

A comment is in order concerning the similarity approach in view of the latest findings 
in the lattice calculations that glueballs are unexpectedly heavy. |Ll| The hypothesis is that 
glueballs are heavy because the effective potentials binding two gluons are much stronger 
(perhaps by a factor of four) than the potentials that bind quarks. The reason is that gluons 
are in the octet representation of SU(3), which is analogous to being a doubly charged object 
in SU(1). Thus, if one compares a quark bound state to hydrogen, a glueball bound state is 
comparable to helium, but with one doubly- charged electron rather than two singly charged 
electrons. Obviously the helium and doubly charged electron are far more strongly bound 
than hydrogen is. 

What this means is that as the coupling constant increases with decreasing hamiltonian 
width, the coupling of gluons will always be four times larger than the coupling for quarks, 
which suggests the gluon coupling would be expected to become strong enough to create a 
bound state of gluons while still too small to bind quarks. The strong gluon binding would 
naturally lead to a high effective mass for gluons, making them unlikely to be present in quark 
bound states. This would explain why we see no evidence of explicit constituent bound states 
involving gluons as well as quarks, yet allow the quark-quark potential to have strong higher 
order corrections. A major question for this picture would be to understand sum rules 
for deep inelastic scattering which have been interpreted as implying a large contribution 
from constituent gluons inside the proton. It might be that development of the similarity 
transformation formalism would show that sum rule data refers more to "current" quarks 
and gluons, before the similarity transformation is applied, rather than to constituent quarks 
and gluons after the transformation, where the gluon contribution could be very small by 
these arguments. 

Our study shows that the convergence of perturbation theory is in jeopardy beyond the 
second order expansion for large values of the running coupling constant. Ref. QXJ] suggested 
that one can use the phenomenological success of the constituent quark model to improve 
convergence when solving QCD in the light-front hamiltonian approach. The improvement 
is expected to result from using a constituent quark hamiltonian as a first approximation to 
the small width effective hamiltonian of QCD. The chance exists, that such a constituent 
quark hamiltonian is not much different from the theoretical one in QCD. Therefore, the 
distance to the true effective hamiltonian may be calculable in perturbation theory. We wish 
to add a comment on how this suggestion could be checked numerically in the model. 

The analogous step could be done in the present model in the following way. For certain Ao 
close to the bound state formation scale, order a few GeV, the running coupling, g = g(\o), 
has a value comparable to 0.5. Let us denote the exact value of go by g s (s is chosen for 
"strong"). The perturbative expansions towards smaller widths order 1 GeV in terms of 
g(X) are hard to continue because of the large distance from the small coupling g in the 
initial hamiltonian. In these circumstances, we can add and subtract a suitable term in 
the hamiltonian TC(\o), say Hcqm- This step changes nothing. But we can multiply the 
subtracted term by the ratio go/g s . For g = g s nothing is changed in the theory. On 
the other hand, if we replace g$ by a small number, the subtracted term together with the 
original interaction can be treated as a small perturbation. In fact, if \xqqm represents the 
bulk of the effective hamiltonian then the difference between the effective hamiltonian and 
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Hcqm will not lead to large corrections even if g is rised to g s . Thus, perturbation theory 
in terms of g could be continued towards the smaller widths and the coefficients could be 
kept small. One can think of the hamiltonian [Icqm as a matrix which has only one matrix 
element different from 0, right in the place on the diagonal where the bound state eigenvalue 
appears. Other forms are also possible. The key example is provided by the approximations 
shown in Fig. 2. The general feature of the example is that it amounts to the insertion of 
the similarity factors and adjusting the couplings. 

Another opportunity for improvement in the numerical accuracy is related to the ir- 
relevant operators. Analysis of the renormalization group equations for low energy matrix 
elements suggests that for small widths A the corrections on the right-hand side in Eq. (2.13) 
should include the term h(X) with matrix elements of the form 

~h(X) 



hmn(X) =^(E m + E n )^E m E n exp [-[E m - E n ] 2 /X 2 \ . (4.1) 

A new coupling constant h ~ g 2 is introduced. The new term behaves as A -1 for large A 
and it disappears for A — > oo. It remains to be verified how much this term can improve the 
accuracy of the perturbative evaluation of narrow effective hamiltonians and the resulting 
bound state eigenvalues. The two couplings g and h are related and should be considered 
coherently. |I2 | 



Too little is known yet about equations of the type (2.14) to say if they can help in 
accelerating convergence of the perturbative expansion. 

We should mention that when b is reduced towards 1 the number of states per unit of the 
energy grows and the resulting matrices would have to be much larger than in the example 
we described. In fact, an interesting problem arises when one considers more than one state 
at each scale. Namely, when we have only one state per scale then no degenerate or near 
neighbour interactions arise. If more states appear at each scale and they couple to the near 
neighbours, the couplings are not reduced by the similarity transformation. The interesting 
question is how large can be the effective energy range of such near neighbour interactions. 
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Figure Captions 

Fig. 1. The approximate running coupling g a (X) from Eq. (2.13) and the exact running 
coupling g(X) from Eq. (3.1) plotted as functions of the effective hamiltonian width A. The 
exact reduced matrix element \x from Eq. (3.2) is also plotted to show the width range where 
the bound state eigenvalue appears on the diagonal. 

Fig. 2. The bound state eigenvalue of approximate effective hamiltonians in ratio to the 
exact value —1 GeV. Curve E results from Eq. (2.13) with corrections set equal 0. Curve 

V is obtained in the same way as E except for H eigenvalues E in the Gaussian similarity 
factor replaced by the diagonal matrix elements of the approximate effective hamiltonian, 
i.e. T> — (1 — g)E, (see the text). The intermediate curve marked g is obtained by using 

V — (1 — g)E in the exponent, i.e. the initial coupling instead of the running one. 

Fig. 3. The ratio of the effective hamiltonian bound state eigenvalue to the exact result, 
for the hamiltonians of width A calculated using expansion in powers of the initial coupling 
constant g up to 1, 2, 3 and 4, respectively. It is visible, that the expansion is not useful in 
the calculation of the bound state eigenvalue. 

Fig. 4. The accuracy of the bound state eigenvalues obtained from effective hamiltonians 
whose renormalization group flow with the width A is calculated expanding in powers of the 
effective coupling constant g{\o) and including terms order 1, g(Ao) and g 2 (Ao). The accuracy 
is given as ratio of the bound state eigenvalue obtained by diagonalization of the effective 
hamiltonian of width A to the exact value, —1 GeV. The curves correspond to the indicated 
values of Ao (in units of GeV). The result of expansion in the initial coupling g is denoted 
by oo. The arrows show points where A = Ao- 

Fig. 5. The accuracy of the perturbatively calculated effective hamiltonians as measured 
by the matrix element jl(X) from Eq. (3.2) in ratio to its exact value. The numbers label 
results obtained by expanding in powers of the running coupling constant g(\) including 
powers up to 1, 2, 3 and 4, respectively. The ratio would be indistinguishable from 1 for 
widths larger than 64 GeV if it were plotted in the scale of this figure. 

Fig. 6. Accuracy of the bound state eigenvalue obtained from diagonalization of the 
effective hamiltonians 7i(A) expanded in powers of the running coupling g(\) including 
powers up to 1, 2, 3 and 4, respectively. The curve 1 matches curve E in Fig. 2 away 
from the left boundary region. The curve order 2 shows 10% accuracy down to the width 
1 GeV. The lack of convergence for higher order curves at small widths is discussed in the 
text. 

Table 1. Ratio of the bound state eigenvalue of the small window hamiltonian with 
indices limited by fh and n, to the eigenvalue of the whole effective hamiltonian at A = 1 
GeV calculated using expansion up to second power in the running coupling g(l GeV). 0.993 
corresponds to the absolute accuracy of the bound state eigenvalue equal 12% and 0.908 to 
19% (see the text). 
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Fig. 1. The approximate running coupling g a (X) from Eq. (2.13) and the exact 
running coupling g(X) from Eq. (3.1) plotted as functions of the effective hamiltonian 
width A. The exact reduced matrix element fx from Eq. (3.2) is also plotted to show the 
width range where the bound state eigenvalue appears on the diagonal. 
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Fig. 2. The bound state eigenvalue of approximate effective hamiltonians in ratio to 
the exact value — f GeV. Curve E results from Eq. (2.13) with corrections set equal 0. 
Curve T> is obtained in the same way as E except for H eigenvalues E in the Gaussian 
similarity factor replaced by the diagonal matrix elements of the approximate effective 
hamiltonian, i.e. T> = (1 — g)E, (see the text). The intermediate curve marked g is 
obtained by using V = (1 — g)E in the exponent, i.e. the initial coupling instead of the 
running one. 
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Fig. 3. The ratio of the effective hamiltonian bound state eigenvalue to the exact 
result, for the hamiltonians of width A calculated using expansion in powers of the initial 
coupling constant g up to 1, 2, 3 and 4, respectively. It is visible, that the expansion is 
not useful in the calculation of the bound state eigenvalue. 
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Fig. 4. The accuracy of the bound state eigenvalues obtained from effective hamil- 
tonians whose renormalization group flow with the width A is calculated expanding in 
powers of the effective coupling constant g(A ) and including terms order 1, g(A ) and 
<? 2 (A ). The accuracy is given as ratio of the bound state eigenvalue obtained by diago- 
nalization of the effective hamiltonian of width A to the exact value, —1 GeV. The curves 
correspond to the indicated values of A (in units of GeV) . The result of expansion in the 
initial coupling g is denoted by oo. The arrows show points where A = Aq. 
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Fig. 5. The accuracy of the perturbatively calculated effective hamiltonians as mea- 
sured by the matrix element jl(X) from Eq. (3.2) in ratio to its exact value. The numbers 
label results obtained by expanding in powers of the running coupling constant g(X) in- 
cluding powers up to 1, 2, 3 and 4, respectively. The ratio would be indistinguishable 
from 1 for widths larger than 64 GeV if it were plotted in the scale of this figure. 
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Fig. 6. Accuracy of the bound state eigenvalue obtained from diagonalization of the 
effective hamiltonians T~C(\) expanded in powers of the running coupling g(X) including 
powers up to 1, 2, 3 and 4, respectively. The curve 1 matches curve E in Fig. 2 away 
from the left boundary region. The curve order 2 shows 10% accuracy down to the width 
1 GeV. The lack of convergence for higher order curves at small widths is discussed in the 
text. 
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Table 1. Ratio of the bound state eigenvalue of the small window hamiltonian with indices 
limited by fh and n, to the eigenvalue of the whole effective hamiltonian at A = 1 GeV 
calculated using expansion up to second power in the running coupling ^(lGeV). 0.993 
corresponds to the absolute accuracy of the bound state eigenvalue equal 12% and 0.908 to 
19% (see the text). 
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